Theory of dark-state polariton collapses and revivals 
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We investigate the dynamics of dark-state polaritons in an atomic ensemble with ground-state 
degeneracy. A signal light pulse may be stored and retrieved from the atomic sample by adiabatic 
variation of the amplitude of a control field. During the storage process, a magnetic field causes 
rotation of the atomic hyperfine coherences, leading to collapses and revivals of the dark-state po- 
lariton number. These collapses and revivals should be observable in measurements of the retrieved 
signal field, as a function of storage time and magnetic field orientation. 
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A quantum memory element consisting of an ensemble 
of atoms, with efficient coupling to a signal light field, 
represents a node in several quantum network architec- 
tures 

nil mi- A dark-state polariton (DSP) is a col- 
lective excitation, with light field and atomic spin wave 
parts, in which the relative size of the light and matter 
contributions can be varied by changing the amplitude 
of a control laser field 1]. In connection with atomic 
memories, DSPs offer the possibility for efficient trans- 
fer of information between a light carrier and an atomic 
medium, with programmable storage of the excitation in 
the atomic spin coherence. The storage and subsequent 
retrieval of the signal field component of the DSP can be 
achieved by the extinction and subsequent reactivation of 
the control field after a given storage time. Experimental 
demonstrations of "stopped-light" can be understood in 
terms of the concept of DSP in just this way 0,0,01- 

In a recent work the storage and retrieval of single pho- 
tons using an atomic ensemble based quantum memory 
was reported, with a storage time conjectured to be lim- 
ited by inhomogeneous broadening in the ambient mag- 
netic field II]. During the storage, the DSP consists en- 
tirely of atomic spin-wave, and in order to understand 
its dynamics in a magnetic field it is necessary to prop- 
erly account for the atomic level degeneracy and the sig- 
nal and control field polarizations. In particular for al- 
kali atoms, which have non-zero nuclear spin, the elec- 
tronic levels have hyperfine structure. In this case we 
must define a more general form of DSP field operator 
than that of a simple lambda configuration, in which the 
atomic spin-wave part corresponds to a particular super- 
position of hyperfine coherences of the ground electronic 
level. These coherences are, in turn, intimately related 
to the phenomenon of electromagnetically-induced trans- 
parency (EIT) We shall see that in a magnetic 
field the temporal evolution of the DSP reveals a series 
of collapses and revivals due to the evolution of its spin- 
wave component during the storage phase of the process. 
We predict that the collapses and revivals should be di- 
rectly observable in measurements of the retrieved signal 
field as a function of storage time. 

We develop the theory of EIT in a degenerate atomic 
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FIG. 1: (Color online) On the left, a diagram shows an atomic 
ensemble interacting with copropagating signal and control 
fields. The signal (helicity a), resonant on the <; <-> e tran- 
sition, is stored and subsequently retrieved by variation of a 
control field (helicity /3), which resonantly couples levels g' 
and e. A constant magnetic field B, oriented at an angle 9 
from the propagation axis, rotates the atomic coherences dur- 
ing the storage. For each state \g, m) in level g, there is an 
associated A configuration, as shown on the right. The signal 
connects the states \g,m) and |e,m -I- a), while the control 
field drives transitions between \e,m + a) and \g' ,m + a — P) . 



medium with ground levels g and g' and excited level e, 
which have energies hwg = 0, hujgi and hiUe, respectively 
(Fig. The Zeeman states of level g are written \g, m), 
where —Fg < m < Fg] similar definitions hold for the 
other levels. All N atoms are assumed to be initially 
prepared in level g without polarization, i.e., the den- 
sity matrix of atom = '^^p\g,'m) ^{g,m\ where 
we write p — l/{2Fg + 1). The density of atoms in the 
sample is assumed to be constant. The atoms experience 
a uniform magnetic field B oriented at an angle 9 with 
respect to the light propagation z axis. The magnetic 
field-atom interaction Vb = A*bB • J2s=g g' e 9s^s, where 
Fs is the projection of the atomic angular momentum op- 
erator onto level s and gs is the corresponding Lande g 
factor. The magnetic field induces a Larmor spin preces- 
sion which is primarily important in the storage phase, 
when the signal field amplitude is zero. We note that in 
prior work collapses and revivals of single-atom Zeeman 
coherences have been observed . 

We proceed by gen eralizing the treatment of Fleis- 
chhauer and Lukin to include the degenerate atomic 
level scheme and the presence of a magnetic field. We as- 
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sume the number of photons contained in the signal pulse 
is much less than the number of atoms in the sample. 
For an initially unpolarized sample in level g, the ground 
state populations and Zeeman coherences, as opposed to 
hyperfine coherences, are unaffected by the signal, con- 
trol and magnetic fields. Our treatment can be extended 
to an initially spin-polarized atomic sample, as we will re- 
port in a separate publication. The signal field, which we 
assume propagates in the positive z direction, and atomic 
coherence operators satisfy the quantum Langevin equa- 
tions 
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where the slowly varying electric field of helicity a = ± 1 
is given by ^aiz,t) = i X]fc a/c,a exp(i((j2; + ujet)), q = 
k — {uje/c) and the collective atomic coherence is de- 
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the /x'th atom hyperfine coherence operator, and = 
Ndz/L is the number of atoms between z and z + dz; L 
is the length of the sample. The control field is assumed 
to have the circular polarization f3 — ±1 and real Rabi 
frequency ^{t), which is a specified function of time, and 
K is the coupling constant for the probe transition. We 
adopt the shorthand for the Clebsch-Gordan coefficients 

Cma = C'ra^a m+a ^^'^ ^'m[3 = p m+(}'^ USeful 

to define Rma{f3) = CmQ/C^_,_„_^_^. The decay rate of 

level e is denoted by and ^, is a corresponding 
quantum noise operator. The coupling of the atoms to 
the uniform magnetic field can be taken into account by 
the addition of appropriate commutators with the inter- 
action Vb in the atomic equations. 

We first establish some standard features of EIT with 
our model. The propagation of a classical (coherent) sig- 
nal through the medium is found by dropping the quan- 
tum noise operator, and replacing the field and coherence 
operators with their respective expectation values. For a 
constant amplitude control field, the linear susceptibility 
for the signal field of angular frequency lu is found to be 
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where A = a; — ojg is the detuning of the signal from 
atomic resonance, and Xraa = C'ma/ \/J2m' ^m'a- The 
dimensionless quantity da is the optical thickness, which 
is defined such that exp(— da) is the on-resonance inten- 
sity transmittance in the absence of a control field, and 
can be expressed as 
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where 77 is the fraction of atoms in excited level e that 
spontaneously decay into ground level g, and A is the 
cross sectional area of the ensemble. When a control 
field is present, an EIT window exists provided that the 
Clebsch-Gordan coefficients C^_|_^_^ ^ do not vanish for 
any —Fg<m<Fg for which Cma 7^ 0. If, however, 
C!^_^_a_p ^ = 0, and Cma 7^ 0, it means that there is an 
excited state \e,m + a) not coupled by the control field to 
a state in the ground level g', i.e., there is an unconnected 
lambda configuration. The subset of atoms initially in 
the state \g, m) would absorb the signal field and spon- 
taneously emit radiation as if there were no control field 
present. In order for EIT to exist, one must make a judi- 
cious choice of atomic levels and signal and control field 
polarizations. 

Assuming a choice of polarizations that supports EIT, 
we are able to generalize the adiabatic treatment of Ref. 
1] to Eq. ^ to derive the DSP operator for helicity a, 
with control field polarization (3 
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As in Ref. Q , this operator obeys the simple propagation 
equation {d/dt + Vgd/dz)^a{z, t) — with the reduced 
group velocity Vg = cn^/{n^ + Np\K\^J2m\R^^W\'') 
which can be adiabatically controlled by time dependent 
variation of n(t). From the definition of we see that 
as f2 goes to zero, the wave excitation stops propagating 
and transforms into a particular linear combination of hy- 
perfine coherences ~ Em -Rm«(/3)Qg'1+a-/3(^' This 
nontrivial result arises from the treatment of the full de- 
generacy of the atomic ensemble; only this combination 
of hyperfine coherences is adiabatically transformed into 
the signal field via the control field retrieval process. Or- 
thogonal combinations of hyperfine coherences couple to 
optical coherences in the presence of the control field and 
result in excited state spontaneous emission; we will refer 
to these as the bright-state polariton (BSP) component. 
It is also possible that some population of atoms remains 
trapped in the ground states, and is unaffected by the 
control field. 

' In order to demonstrate the importance of the dark 
state polariton in the signal storage and retrieval process 
in a magnetic field, we numerically solve Eqs. (^3) for 



a coherent signal field. We thus calculate the expecta- 
tion values of the spin wave coherences {Q^^j^i^_p{z, t)) 
and the signal field, allowing us to determine the DSP 
and BSP components. This is accomplished by defin- 
ing a vector space of 2Fg + 2 dimensions, with or- 
thonormal basis vectors e^, each corresponding to a 
hyperfine coherence and e$ corresponding 

to the signal field. We define the coherence "vector" 

V ^ (^o)e*-KE™y^(Qyi+„-^(^,t))e™. We 
note that this is not normalized since its magnitude is 
dependent on both the time dependent signal and con- 
trol fields. Associated with the DSP we define a vector 
uq, = rie$ + \/NpK J2m Rma{f3)em and the correspond- 
ing unit vector e^i = u$/||u$||. We then determine the 
DSP component pD = |e^ • vp, and BSP component 
PS = ||v - e^e* • v|p. 

As an example of signal storage and retrieval we con- 
sider an atomic sample of *^Rb, in which the control field 
and signal field polarizations are chosen equal a = f3 = 1, 
and the optical thickness d = 8, Fig. [3 The atomic levels 
g,g' and e correspond to the 5Si/2 F — 2,3 and 5Pi/2 
F = 3 levels of the Di line, respectively. The sponta- 
neous decay rate re/(27r) = 5.98 MHz. The incident sig- 
nal field has a Gaussian envelope of full width half max- 
imum 120 ns, and the peak enters the 3 mm long sample 
at t = —60 ns. The control field has a constant Rabi 
frequency — l.SFe until it is smoothly turned off at 
t = over a period of 20 ns, when a fraction of the signal 
field is converted into hyperfine coherences of the atomic 
spin wave. The excitation is stored from < t < 2 /is in 
the presence of the magnetic field, before the control field 
is reactivated, and the signal field retrieved. In Fig. El 
panels (A) and (B), the magnetic field is chosen so that 
the storage time corresponds to a quarter of a Larmor pe- 
riod Tl = 2Trh/\ggfiB^\, while in panels (C) and (D), the 
storage time is Ti/2. In panel (A) pn grows as the signal 
pulse arrives at the point of observation, and reaches a 
peak when the control field is switched off. It then decays 
during the storage phase, due to Larmor precession of the 
hyperfine coherences in the applied magnetic field, which 
causes the corresponding growth of pb ■ When the control 
field is reactivated, pb decays rapidly due excited level 
coupling and subsequent spontaneous emission, though 
Pd remains finite as the spin wave coherence of the DSP 
is converted into the forward propagating signal field; the 
retrieved signal field intensity is illustrated in panel(B). 
In panel (C), where the storage time is Tl/2, pd un- 
dergoes a complete revival. The energy of the retrieved 
signal field shown in panel (D) is therefore much larger, 
by a factor of 5.73, than that in panel (B). This is in 
good agreement with the DSP theory for retrieval effi- 
ciency discussed later, which predicts that the retrieved 
signal energy of panel (D) should be 5.53 times that in 
panel (B). These results demonstrate the importance of 
the adiabatic concept of DSP for a realistic experimental 
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FIG. 2: (Color online) Numerical results illustrate the stor- 
age and retrieval of a signal pulse from an atomic ensemble as 
described in the text. We show results for two values of the 
magnetic field oriented along the z axis. Panels (A) and (B) 
correspond to a magnetic field B — 0.267 G, so that the sig- 
nal is stored for Tl /4, where Tl — 8 ^is is the Larmor period. 
Panels (C) and (D) show results for B — 0.535 G, correspond- 
ing to a signal storage time of Tl/2, where Tl = 4/^s. The 
signal field intensity transmittance 7(t)/7o (Solid line) and 
control field Rabi frequency (dot-dashed line), displayed in 
arbitrary units, are shown in panels (B) and (D). Panels (A) 
and (C) display scaled dark state po (solid line) and bright 
state Pb (dashed line) polariton components, as explained in 
the text. In panel (B) the ratio of retrieved to input signal 
pulse energy is 4.38% while in (D) the ratio is 25.09%. 

scenario. The retrieved signal field directly reflects the 
DSP dynamics in the magnetic field. 

We can predict the retrieval efficiency of a stored signal 
pulse by tracking the population of the DSP as it evolves 
under the infiuence of the magnetic field. To compute 
the polariton population, it is convenient to consider the 
Fourier components of the DSP. We express the polariton 
annihilation operator for the mode of wavenumber q as 

^Jn^ + Np\K\'j:mR^c{f3y 

(5) 

where S*^,™, = (l/\/iV^) '^Jm.g'm' exp(-i(<72;^ + 
{ujg — ujgi){t — z^/c))) is a collective spin wave 
annihilation operator of wavenumber q. These op- 
erators obey quasi-bosonic commutation relations 

the DSP operators therefore also obey [^'a(g), ^'|,/ (g')] — 
Sqq'Saa' + 0{1/N). During the storage, the evolution 
of the spin wave operators is given by S^l^^_p{q,t) — 

where 2?,|^|„/(i) = {s,m\exp{-igsflB ■ Ft)\s,m') is 
the matrix element of the rotation operator for states 
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FIG. 3: (Color online) The DSP population fraction fct{ts) 
calculated for orientations of the magnetic field < 9 < n /2 
over one Larmor period. These results illustrate collapses 
and revivals whose features are dependent on 9. The atomic 
configuration and field polarizations a and /? are described in 
the text. 



in hyperfinc level s, and Ob = /ibB/S. Using the 
bosonic commutation relations for the spin wave op- 
erators, we can calculate the number of polaritons 
{Na{ts)) = (Eg *l(9,is)*a(9,is)> as a function of stor- 
age time ts for an arbitrary DSP quantum state created 
in the storage process. In the limit of infinite control 
field amplitude, this converts into the total number of 
photons in the retrieved signal field J2k '^k a^k,a- We 
therefore derive an expression for the signal retrieval 
efficiency as the fraction fa{ts) = {Na{ts)) / {Na{0))'- 
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In Fig. 121 we display the fa{ts) as a function of tg for 
a variety of magnetic field orientations. We again con- 
sider an ensemble of ^^Rb atoms with the same choice 
of atomic configuration and field polarizations discussed 
earlier. For tg <C II, we observe a collapse in the po- 
lariton population, yielding an approximate retrieval ef- 
ficiency of fa{ts) ~ exp(— 7y^(rJLis)^/2), where the col- 
lapse rate r/a depends on the angle, 9, between the mag- 
netic field and the propagation axis. For ^ = 0, we find 



thus the signal retrieval efficiency should equal the zero 
storage time value. Depending on the orientation of the 
magnetic field, we observe also a partial revival at half 
the Larmor period. For a magnetic field oriented along 
the z axis, the system dynamics are relatively simple. 
Each hyperfine coherence merely picks up a 

phase factor that oscillates at m + {a — P)/2 times twice 
the Larmor frequency, thus returning the system to its 
initial state at half the Larmor period. In this case, the 
partial revival is actually a full revival. On the other 
hand, for = 7r/2, a rotation through half the Larmor pe- 



riod causes the coherence transformation S^, 



As a re- 



up to an overall phase factor, 
suit, for the choice of equal field polarizations (a = 
the retrieval efficiency at half the Larmor period simpli- 
fies to (I],„i?ma(Q:)i?-m,a(a)/I]„ |^?mQ(a)P)^ result- 
ing in a partial revival. For other orientations of the 
magnetic field, particularly for 9 — 7r/4, the revival at 
half the Larmor period is suppressed. This refiects the 
more complicated dynamics of the individual spin coher- 
ences Sg, ™j, each of which transform into a superposition 
of all {2Fg + l){2Fg' + 1) spin coherences, with complex 
time dependent coefficients governed by the rotation ma- 
trices. Stated physically, there is a strong destructive 
interference between the various spin coherences when 
9 Kin/A. 

We have developed a theory of the DSP as a mech- 
anism to store and retrieve light pulses in a degenerate 
unpolarized atomic medium. The role of the DSP and 
its connection to storage retrieval efficiency have been 
verified by full numerical solutions of the propagation 
equations for a classical incident signal field. In the pres- 
ence of a magnetic field, we have demonstrated that the 
DSP population undergoes collapses and revivals during 
the pulse storage time. We predict that this polariton 
dynamics is directly reflected in the signal pulse retrieval 
efficiency. Our results may find applications in quantum 
communication and computation approaches that utilize 
quantum memories [3, UM UM ■ 

This work was supported by NASA, Office of Naval 
Research, National Science Foundation, Research Cor- 
poration, Alfred P. Sloan Foundation, and CuUen-Peck 
Chair. 

Note added: While this manuscript was beiiig pre- 
pared, the collapses and revivals were observed [l9|, in 
excellent agreement with the predictions of this work. 
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With the approximation gg ~ —gg' , valid for ground level 
alkalis, it is clear that the system undergoes a revival 
to the initial state after a complete Larmor period, and 
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